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THE  EFFECT  OF  THICKNESS  ON  PITCHING  AIRFOILS  AT  SUPERSONIC  SPEEDS 


ABSTRACT 

The  effect  of  thickness  on  the  lift  and  pitching  moment  on  two 
dimensional  pitching  airfoils  at  supersonic  speeds  is  investigated. 

The  airfoils  considered  have  arbitrary  symmetrical  cross  sections,  and 
the  flow  is  supersonic  throughout  the  flew  field* 

The  analysis  is  based  on  a second  order  theory  similar  to  the  second 
order  theory  introduced  by  Busemann  and  extended  by  Van  Dyke.  The  lifting 
pressure  due  to  steady  pitching  is  found  and  this  is  used  to  calculate  the 
lift  due  to  pitching,  CT  , and  the  moment  due  to  pitching,  C . 
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INTRODUCTION 


The  development  of  the  linearized  theory  of  supersonic  flow  ..as 
permitted  a first  order  evaluation  of  a number  of  stability  derivatives* 
Second  order  theories  similar  to  the  one  introduced  by  Busemann  (ref.l) 
and  extended  by  Van  Dyke  (ref. 2 and  3)  offer  possibilities  of  obtaining 
second  order  evaluations  of  certain  stability  derivatives*  such  as  lift 
and  moment  due  to  steady  Ditching,  and  lift  and  momer'.  due  to  constant 
rate  of  change  of  angle  of  attack.  The  determination  of  the  dairying  in 
roll  (in  ref.  I4)  for  certain  airfoili  is  an  example  of  the  use  of  a second 
order  theory  to  obtain  stability  derivatives. 

In  this  paper  a second  order  theory  is  developed  for  two  dimensional 
pitching  airfoils  at  supersonic  speeds.  This  theory  yields  an  expression 
for  the  lifting  pressure  due  to  steady  pitching  which  enables  the  stability 
derivatives  CT  (lift  due  td  pitching)  and  C (moment  due  to  pitching)  to 

L |Q 

be  calculated.  The  airfoils  considered  here  ^have  arbitrary  symmetrical 
cross  sections;  however,  the  analysis  can  easily  be  extended  to  include  air- 
foils with^un symmetrical  cross  sections. 

The  partial  differential  equation  for  airfoils  with  a steady  pitching 
velocity  is  expressible  in  a form  independent  of  time.  Work  by  Milton  D. 

Van  Dyke  (ref.  2 and  3)  indicates  that  second  order  solutions  of  the 
* partial  differential  equation  of  steady  supersonic  flow  can  be  obtained  by 
iterative .methods.  The  equation  considered  here  is  quite  similar  to, 
though  not  the  same  as,  the  equation  of  steady  supersonic  flaw.  We  *all 
assume  that  the  second  order  solution  can  be  obtained  by  the  use  of 
iterative  methods. 

It  will  also  be  assumed  that  the  characteristics  are  the  saJoe  for  the 
first  and  the  second  order  solutions  (this  assumption  was  made  ih  ref.  2 
and  3).  For  steady  plane  flow  the  second  order  solution  (ref.  2)  found  by 
using  an  iterative  method  based  on  the  preceding  assumption  yields  the 
correct  second  order  pressure  of  the  Busemann  second  order  theory.  Un- 
fortunately, no  such  justification  of  this  assumption  is  known  to  the 
authors  for  the  flow  associated  with  pitching  airfoils*. 

ANALYSIS 


The  Partial  Differential  Equation:  The  partial  differential  equation 

to  be  used  in  the  following  analysis  is  a special  case  of  the  two  dimensional 
time  dependent  equation  for  the  potential  function  of  a non-viscous  compressible 
fluid.  This  equation  is  (see  ref.  U ) 
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This  equation  is  associated  with  axes  fixed  in  space.  Far  the  present 
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problem  it  Is  convenient  to  express  the  differential  equation  in  terms 
of  a set  of  axes  which  are  fixed  to  the  airfoil.  Since  an  airfoil  with 
a constant  rate  of  pitch  can  be  considered  as  flying  in  a circular  flight 
path  with  a constant  speed,  the  axes  fixed  to  the  airfoil  must  rotate  yith 
respect  to  the  stationary  axes  with  a constant  angular  velocity  q.  The 
relations  between  the  two  sets  of  coordinates,  which  are  illustrated  in 
fig.  1,  are 


x - xQ  cos  qt  ♦ aQ  sin  qt  (a) 

z - V/q  - xq  sin  qt  «►  cos  qt  (b)  (2) 

The  flow  over  an  airfoil  with  a constant  rate  of  pitch  is  steady,  thus 
in  the  new  coordinate  system  the  partial  derivatives  with  respect  to  time 
will  be  zero  in  the  present  problem.  In  the  axes  attached  to  the  airfoil 
equation  (1)  i3  to  the  second  order 

-2  *xx  - [*2<r  - ♦ 2M2q*#  M *M2  - 

ZK2qz  $ ^ 2M2  ] / T (3) 


The  Pressure  Relation:  The  tine  dependent  pressure  relation  is  (from 

ref.  1) 

% * - [2(#Io  A)  /']  -(#,2  * 

In  the  axes  attached  to  the  airfoil  the  preceding  equation  becomes 

cp  A2)"^z2A2j+^2  * 

Eqs.  (3)  and  (U)  will  be  used  in  the  following  analysis. 


Solution  by  Iteration:  It  will  be  assumed  that  eq.  (3)  can  be  solved 
by  an  iteration  procedure.  The  first  order  solution  is  taken  as  the  first 
approximation  to  the  solution  of  eq.  (3).  The  first  order  differential 
equation  i3  obtained  by  neglecting  the  second  order  tfcrms  in  eq.  (3),  arxi 
it  is  given  by 
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It  is  assumed  that  the  second  approximation  can  be  found  by 
substituting  the  first  order  solution  into  the  right  side  of  eq. 

(3)  and  solving  the  resulting  nonhomogeneous  equation,  which  from 
eq.  (3)  and  (5)  is 

-P2  *xx+  Kz  " [m2(Y  "^x^xx*  *zz>  ♦ 2 q x <*xz  + 

T^q  lfz  - 2M2  q z ♦ 2H2  ^ ^ ♦ 

2M2  0y  ♦ 2M2  jf,  ?„]  A (6) 

The  solution  of  eq.  (6)  will  be  referred  to  as  the  second  order 
solution. 

TH nH nation  of  Terms:  It  is  helpful  to  investigate  the  type 

of  solution  obtained  from  eq.  (6).  The  first  order  equation  will  be 
of  the  form 

0 - o g1  (x,z)  ♦ q g2(x,z)  ♦ £g3(x,z) 

where  £ is  a thickness  parameter  and  where  a,  q,  and  6 are  small  compared 
to  unity.  It  follows  from  the  preceding  expression  and  eq.  (6)  that  the 
second  order  solution  will  be  of  the  form 

irm  a 2 hjCxjz)  + q2h2(x,z)  ♦ * aqhj^x,*)* 

a * h^(x,z)  ♦ q£h^(x,z) 

The  terms  a2h^(x,z),  a^^XjZ)  ^ and  a e h^(x,z)  are 

independent  of  the  rate  of  pitch;  thus  they  do  not  contribute  to  the 
lifting  pressure  due  to  the  pitching  motion  of  the  airfoil  • Since  this 
paper  is  concerned  only  with  the  lifting  pressure  due  to  the  pitching, 
all  terms  of  the  forms  a h(x,z)  , *zh(x,s)  , and  a«h(x,z) 

can  be  neglected  in  the  ensuing  analysis. 

n 

The  following  argument  shows  that  q h2  (x,  z)  ant)  aqhj^Xjz) 

do  not  contribute  to  the  lifting  pressure.  The  thickness  parameter,  £ , 
is  not  present;  therefore,  the  airfoil  can  be  considered  a flat  plate 
insofar  as  these  terms  are  concerned.  If  the  potential  of  the  flow  on 
the  upper  surface  of  a flat  plate  is  expressed  as 

£(<*>  q)  ■ ♦ qg2  + a2hi  + cc<i*V'  q.^2 

the  potential  on  the  lower  surface  is  given  by 

#(-  a,  - q)  - - ogx  - qg2  + ♦ aqh^  ♦ q2^ 

the  potential  difference  is 

• 2 a g1  * 2 <s2 

* The  terms  e^h^x,  z),  and  a(L(z,  z)  are  associated  with  steady  two 
dimensional  supersonic  flow  ana  can  be  found  by  use  of  eq.  (U6)  in 
ref.  (2) 
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Since  fear  the  flat  plate  the  pressure  difference  between  the  upper  and 
lower  surface  can  £e  found  directly  Croa.  the  potential  difference  the 
terns  a q h^  and  q*  h2  do  not  contribute  to  the  lifting  pressure* 

The  remaining  tern  qah^  (x  ,*)  will  be  found  by  use  of  eq*  (6), 
neglecting  all  expressions  multiplied  a2,  a<,  q , a2,  and  aq. 

A further  consideration  of  the  form  of  the  solution  indicates  that 
the  second  order  lifting  pressure  is  linear  in  the  thickness  parameter. 
This  can  be  established  by  considering  a first  order  solution  of  the 
fora 


*)♦  (*,*) 

The  second  order  lifting  pressure  will  be  of  the  form 

ACp  “0*  q *>>  W + ^2  l2 

This  equation  is  linear  in andaL;  thus  the  lifting  pressure  for 
various  known  thickness  distributions  can  be  added  to  obtain  the  lifting 
pressure  for  new  thickness  distributions.  For  the  two  dimensional  airfoil 
this  is  not  of  mud:  value  for  determining  analytical  solutions,  since  we 
shall  determine  the  solution  far  an  arbitrary  thickness  distributioft. 

This  linearity  is  also  true,  however,  for  three  dimensional  airfoils,  and 
for  these  it  should  prove  quite  useful. 

Boundary  Conditions t Physical  considerations  require  that  the  flow  be 
tangent  to  the  surface  of  the  airfoil,  and  that  all  velocity  perturbations 
vanish  upstream  of  the  airfoil.  These  boundary  conditions  may  be  expressed 
mathematically  as: 

0 (x,  z)  - 0 

V (x,  z)  - 0 


upstream  of  the  airfoiA 
leading  edge  Mach  sheet 


and 


oo  .ys  - D 

where  s(x,z)  - 0 is  the  equation  of  the  surface  of  the  airfoil. 

The  equation  of  t he  surface  of  the  airfoil  may  also  be  expressed  as 
z »£f(x) 


Thus 

^s  - - leef/ax  ♦ k 

Since  the  velocity,  co,  can  be  written  as 

co  - i(V  - sq  ♦ ♦ J^x)  ♦ K(qx  ♦ ^ ♦ fz) 

it  follows  that  the  boundary  condition  an  the  body  surface  is 
-(V  -sq  ♦ 0X  +tx)  £fx  ♦ ffz  +£  ♦ qx  - 0 
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Thus  the  boundary  conditions  for  the  first  order  solution  are 
0(x,z)  * 0 upstream  of  the  airfoifs  leading  edge  Mach  sheet,  and 

^.1  * - q 3c  * Vff  f (6) 

•lz-0 

The  boundary  conditions  for  the  second  order  solution  are 
tM  m 0 upstream  of  the  airfoil ' s leading  edge  Mach  sheet,  and 
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zz 


z 
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(9) 


For  the  airfoils  considered  here  the  first  order  velocity  components 
are  discontinuous  across  the  Mach  sheet  from  the  leading  edge. 


To  evaluate  the  effect  on  the  second  order  solution  of  the  dis-  ‘ 
continuities  in  the  first  order  velocity  components,  ve  shall  assume  the 
Mach  sheet  from  the  leading  edge  to  have  a small  thickness.  Within  this 
Mach  sheet  the  first  order  velocity  components  will  be  made  continuous 
so  that  the  discontinuities  through  it  are  replaced  by  continuous 
functions.  This  process  is  illustrated  in  fig.  (2)  for  the  velocity 
component  in  the  free  stream  direction.  The  effect  of  the  -diecontintdtiee 
on  the  r°cond  order  solution  will  be  found  by  obtaining  the  second  order 
solution  within  the  Mach  sheet  and  then  letting  the  thickness  approach 
zero.  This  is  the  same  procedure  followed  in  ref.  (U)  for  rolling  airfoils. 

Only  the  Mach  sheet  above  the  airfoil  will  be  considered  since  the  Mach 
sheets  above  and  below  the  airfoil  are  of  the  same  form.  Since  they  depend 
only  an  the  initial  slope  of  the  airfoil,  the  discontinuities  in  the  first 
order  velocity  components  through  the  leading  edge  Mach  sheet  above  the 
airfoil  will  be  the  same  as  those  through  the  Mach  sheet  above  a flat 
pitching  airfoil  with  a constant  angle  of  attack  of  the  amount  - £, 


The  first  order  potential  function  for  the  flew  over  the  upper  surface 
of  a flat  pitching  airfoil  at  an  angle  of  attack  - € is  given  by 

{*  - [q(x  - pz)2/2  -V6(x-  pz)  - qd2/2  - VS  d]  / p (10) 

where  the  airfoil  is  approximately  in  the  z ■ 0 plane. 

The  discontinuities  in  the  first  order  velocity  components  through 
the  leading  edge  Mach  sheet  above  the  airfoil  are  from  eq.  (10). 

A0X  — q*/P  - V€/p 

A0Z  - q<i  ♦ v« 

It  will  be  assuned  that  the  Mach  sheet  has  a small  thickness  (see  fig.2-b). 
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a.  Leading  Edge  ‘Mach  Sheet 
With  Zero  Thickness 


c.  Discontinuity  In  +,  Across 
Leading  Edge  Mach  Sheet 
Along  The  Line  Z*Z| 


z 


b.  Leading  Edge  Mach  Sheet 
With  Assumed  Thickness 


d.  Plot  Of  4,  Across  Leading 
Edge  Mach  Sheet  With 
Assumed  Thickness  Along 
The  Line  Z*Z( 


Fig  2 -An  illustration  of  removing  the  discontinuity  in 

across  the  leading  edge  Mach  sheet  by  assuming 
the  Mach  sheet  to  have  thickness. 
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The  velocity  components  within  thB  sheet,  will  be  defined  as 
- -(qd  +?«  )(x-p*+d  + f)/  ( f p) 

<fz  * (qd  +V€)  (x-ps+d*/)/  f (11) 

where  f is  the  thickness  of  the  Mach  sheet  in  the  x direction*  Note  that 
within  the  Mach  sheet 

- B2  $ * d - 0- 

h *xx  . *zz 

and  that  the  velocity  components  are  continuous  functions  in  the  neighbor- 
hood of  it* 


From  eq*  (6)  and  (11)  the  differential  equation  within  the  Mach  sheet 
5 s 


-e2 


2 dM1‘(  Y ♦ 1)  ne(x  - p I + d + f)/  (p2  f2)]  * 
2 M2  q«p  z/(p2/)]*-  [M2qc  (x  - p*  ♦ d ♦ f)  / /] 


It  is  convenient  to  express  eq,  (12)  in  terms  of  variables  which  lie 
along  the  characteristics.  Let 

<J-  x - p z 

7]  - x ♦ p z 

In  terms  of  the  new  coordinates  (£,  rj  ) eq,  (12)  becomes 

V -[W*(Y  * 1)  q«  (£♦  d * f)/(U  p1*  f2)]  -[m2**  (£♦  d ♦ (D/Oijfp2)]  - 
[m2  q<  (£*7)  /(!»  p2/)J -[m24«(7-0/(UP1‘  £)] 

This  equation  can  be  integrated  to  yield 

- ^drt^y  ♦ 1)  q«/  (U  PU  t2)] & d ♦ f)  d<*  d^  - 
_M2q €/(U  p2  ^ d ♦ f)  d<*dy-  [M2  qc/(U  p? 

[M2q€/(UpUjT )] d£d  7/ 

The  potential  on  the  downstream  side  of  the  Mach  sheet  can  be  found  by 
evaluating  the  integrals  in  the  preceding  equation  and  taking  the  limit  as 
jT—O. 


O^dCdij- 


The  value  of  the  discontinuity  in  the  second  order  jxrtential  function 
across  ^he  leading  edge  Mach  sheet  is  (from  the  above  equation), 

■ - M2  q€  f * (-d)  ("m2x  ♦ (2  ♦ yM^dl  pz/(2pli) 

- flz  — d L J 


(13) 
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Solution  of  the  Partial  Differential  Equation t The  pert  of  the 
second  order  potential  function  which  yields  the  lifting  preastre 
due  to  the  pitching  moticu  will  now  be  determined.  The  determinations 
of  the  potentials  over  the  upper  and  lower  surfaces  are  similar}  there- 
fore, only  the  flow  over  the  upper  surface  will  be  considered' in  detail. 

The  first  order  solution  is 

0 - [q  (x  - pz)2/(2  P)]  - [v«f(x  - pz)/  p]-  qd2/(2  p)  (li*) 

It  follows  from  eq.  (6)  and  (ill)  that  the  second  order  potential 
function  must  satisfy  the  nonhonogeneous  equation 

- pVxx*  ?zz  ’ {f’(x  - Pz)  - [Ar  ♦ D/P2J|(*  - pz)f"(x-p*)*^ 

f’(x  - pz)]*  2 (px  ♦ z)  f"  (x  - pz)  / p] 

where  the  primes  denote  derivatives  with  respect  to  (x  - pz),  and  the 
a,  aq,  a 6 , and  q^  terms  have  been  neglected^  By  inspection  a partic- 
ular solution  of  eq.  (15)  is  found  to  be 

fx  - «2  {-  [2  -p  2 )/(2  P)]+[m2(y  * i)*(x  -pz)f'(x-pz)/(2p3)]t 

(16) 

+ jj-2  x pz  f’(x  - pz)  ♦ (p2-  1) z2f ' (x-pz )*(p2-l)zf (x-pz)/pj  /(2p2)]- 


The  solution  of  eq.  ( 1 consists  of  a particular  integral, ifr., 
plus  a complementary  function,  |K,  which  must  be  found  so  that  if 
satisfies  the  boundary  conditions  given  by  eq*  (9)  and  (lU).  Hence  ifr^ 
must  satisfy 


- f>2  (f2)xx  * «$>„  - 0 


and  the  following  boundary  conditions: 


•t 

pz  - d x 


pz 


X 


pz  - d 


a*2^z  | z - 0 " ! z 


o -9>V^Z 


/ 


(17) 


- [W(??3)]  (-W1*  T *(pli*2p2-l)xf,(x)^L+p2-2p1+)f(x)](l8) 


By  inspection  the  ccnplejnentary  function  is  found  to  be 


1$ 


The  desired  sojLu^ion  is 

^-*1  *t2  - [q<-/(2pU)J  [mV  Y + 1 - p2)  ? x » f (z  - p s)  - 
M2(M2y  - 2)  p2  t2  f(x  - pz)  + F -M2  pz  ♦ y<*  - P»)  - 

(pU  + 2 p2-  1)  (x  - pz)]  f(x  - pz)  + [-M^  y + 3 p^p2-  ^ f 

Eq'«  (20)  is  the  part  of  the  second  order  potential  function  which 
yields  the  entire  second  order  contribution  to  the  lifting  pressure  due  to 
pitching. 

Lifting  Pressure;  The  lifting  pressure  distribution  can  be  expressed 
as 


C 

P 


upper 

surface 


lower 

surface 


From  eq.  (h)  and  (20)  it  follows  that  for  symmetrical  airfoils  the  lifting 
pressure  distribution  is 

&cp  - - [Uq/(pv)j[x  +^/(2p3)j  |(M^y  ♦ P^-2p2*  l)xf'(x)-H2f(x)|J(21) 

Fig.  (3)  presents  the  lifting  pressure  distribution  on  a ten  per  cent 
thick  wedge  pitching  about  the  c/2  point  for  various  Mach  numbers. 

Fig.  (U)  presents  the  pressure  distribution  on  a ten  per  cent  wedge  at 
Mach  number  1.5  for  various  positions  of  the  axis  of  pitch. 


Fig.  (5)  presents  the  lifting  pressure  distribution  on  a five  per  cent 
thick  airfoil  with  a parabolic  cross  section  pitching  about  the  c/2  point 
for  various  Mach  numbers. 

Fig.  (6)  presents  the  lifting  pressure  distribution  on  a five  per 
cent  thick  airfoil  with  a parabolic  cross  section  at  Mach  number  1.5  for 
various  positions  of  the  axis  of  pitch. 


Stability  Derivatives?  The  effect  of  thickness  on  the  lift  and 
moment,  and  Cm  , due  to  steady  pitching  can  be  found  by  use  of  eq.  (li). 

The  stability  derivatives  and  Cm  can  be  expressed  as 


_q  mq  ~ 

- ^CL/9(qc/(2V)|j  q-*0  - [l/ |qc2/(2V^J_< 

l 


-d 

c-d 


Acp  dx 


x Ac  dx 

P 
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M*l.5 


Percent  Of  Chord 

Fig.  5 The  lifting  pressure  distribution  on  o five  percent  thick 
airfoil  with  parabolic  cross  section  pitching  about  the  c/2 
point  for  Mach  numbers  1.5  and  2.5  . 
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Fig  6 The  lifting  pressure  distribution  on  a five  percent  thick' 
airfoil  with  parabolic  cross  section  at  Mach  number  1.5 
for  various  positions  of  the  axis  of  pitch. 


Fran  the  preceding  relations  and  eq.  (21)  the  C_  and  the  C of  an 
airfoil  are  Lq  ®q 

C.  - Bd/(pc)  - U/p  Y *(p2-l)2Jfl^i/c)/c]  tjg  - 

|’nk  T ♦(P8-1)2  ♦ M2]^  [f(<£)/cZj  d<| 


and 


C,  - - [8/(3p  jj  [l  - 3 d/c  ♦ 3 d2/c2  ] - 

(Uf/p^fyv  ♦(p2-d2j  [d-d/c)2/c]  fra  - 

[W4  Y ♦ 2(p2  l)2  ♦ M 2jJ^  [<lf(4  )/c^|  d<*j 


(23) 


Bq.  (23)  indicates  that  the  effect  of  thickness  on  theX  is  aero  for  an 
airfoil  pitching  about  the  e/2  point  with  a thickness  ^^distribution 
which  is  cHord-Wise  symmetrical  with  respect  to  the  mid-cnord  point* 


Fig.  7 through  10  present  the  variation  of  thA  Cr  for  a ten  f&e  cent 
thick  wedge  and  a five  per  cent  parabolic  airfoil  wit  £4  Mach  number"  and  the 
position  of  the  axis  of  pitch*  Tig*  7 and  9 indicate  that  for  d » 0 the 
thickness  decreases  the  oft  wedge  airfoil  and  increases  the  0^  far 

a parabolic  airfoil*  **  ** 


Fig.  11  through  lii  present  the  variation  of  the  C for  a ten  per  cent 
thick  wedge  and  a five  per  cent  parabolic  airfoil  vit  ft\q/  Mach  number  and 
the  position  of  the  axis  of  pitch*  Fig*  11  and  13  Indicate  that  for  d - O 
the  thickness  decreases  the  far  a wedge  airfoil  and  increases  the  C 

far  a parabolic  airfoil*  ^ *** 

coNCLUnruo  rjwahxs 


The  airfoils  considered  in  this  paper  have  symmetrical  thickness 
distributions*  But  since  the  flow  over  the  upper  and  lower  surfaces  of 
the  airfoils  treated  hdre  are  independent  of  each  other,  the  aerodynamic 
properties  due  to  pitching  of  airfoils  with  un symmetrical  thickness 
distributions  can  easily  be  determined  from  the  results  obtained  here* 

The  limitations  of  the  Busesann  second  order  theory  have  been 
investigated  (see  ref*  5)  Since  the  theory  contained  in  the  present  paper v 
is  closely  associated  with  the  Busemaxm  second  order  theory’,  it  seats  likely 
that  the  results  presented  herein  have  similar  limitations* 


C_  v \\A 

JOHH-JC.  hasten 

NATHAN  GERBER 
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15  2.0  2.5  3.0 

Mach  Number -M 

Fig  7 The  variation  of  the  Ci_q  of  a ten  percent  thick  wedge  with 
Mach  number  for  various  positions  of  the  axis  of  pitch. 
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Fig  8 The  variotion  of  the  Ci_q  of  a ten  percent  thick  wedge 
with  the  location  of  the  axis  of  pitch  for  Mach  numbers 
L5  and  2.5. 


2; 


Fig  9.  The  variation  of  the  Ci_q  of  a five  percent  thick  parabolic 
airfoil  with  Mach  number  tor  various  positions  of  the  axis 
of  pitch. 


Fig  10  The  voriotion  of  the  CLq  of  a five  percent  thick  parabolic 
airfoil  with  the  location  of  the  axis  of  pitch  for  Mach  numbers 
1.5  and  2.5. 


Fig  12  The  variation  of  the  Cmq  of  a ten  percent  thick  wedge 
airfoil  with  the  location  of  the  axis  of  pitch  for  Mach 
numbers  1.5  and  2.5. 
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Fig.  13.  The  variation  of  the  Cmq  of  o five  percent  thick  parabolic 
airfoil  with  Mach  number  for  various  positions  of  the  axis 
of  pitch. 
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FIG.  14.  The  variation  of  the  Cmq  of  a five  percent  thick  parabolic 
airfoil  with  the  location  of  the  axis  of  pitch  for  Mach 
numbers  1.5  and  2.5. 
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